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Spin wave dispersion and damping are investigated in the metallic SDW state of 
different itinerant electron models including a small interlayer hopping. Magnetic 
excitations in iron pnictides are shown to be well understood in terms of physical 
mechanisms characteristic of metallic magnets, such as carrier-induced ferromagnetic 
spin couplings, intra-band particle- hole excitations, and the spin-charge coupling 
mechanism, which is also important in ferromagnetic manganites. 
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I. INTRODUCTION 

Single-crystal neutron scattering studies of iron pnictides have indicated a commensurate 
magnetic ordering of iron moments ordered ferromagnetically in the b direction and antifer- 
romagnetically in the a and c directions.- Inelastic neutron scattering measurements^- in 
AFe2As2 (A = Ca, Ba, Sr) yield well-defined spin-wave excitations up to the zone bound- 
ary on an energy scale ~ 200meV. The realization of a (tt, 0, tt) ordered SDW state has 
opened the possibility of observing phenomena in this class of compounds which are char- 
acteristically associated with both the antiferromagnetic (AF) state such as quantum spin 
fluctuations, hole/electron motion in AF background, spin-fluctuation mediated pairing, 
as well as the metallic ferromagnetic (F) state such as carrier-induced spin interactions, 
correlation-induced spin-charge coupling, and non-quasiparticle states. 

The metallic SDW state in principle also allows for spin wave damping due to decay 
into low-energy intra-band particle-hole excitations, in addition to the usual gapped inter- 
band (Stoner) excitations. Indeed, from the observed high energy behaviour of spin wave 
damping ascribed to particle-hole excitations,- it has been inferred that the full excitation 
spectrum can not be understood in terms of the localized spin models." Weakly damped 
spin waves near the ordering wavevector have been obtained within multiband models from 
the imaginary part of the spin fluctuation propagator .— ~— Spin wave energy renormalization 
and spectral broadening have also been investigated in the metallic AF state.— 

Spin wave dispersion in the (0, tt) and (0, tt, tt) ordered SDW states and the role of doping- 
induced ferromagnetic spin couplings mediated by the exchange of the particle-hole prop- 
agator has been investigated within single-band and two-band^ models with the Hubbard 
interaction^— and with a local spin-fermion exchange coupling^ as in the ferromagnetic 
Kondo lattice model for ferromagnetic manganites. 

The presence of ferromagnetic chains in the SDW state of iron pnictides also suggests the 
possibility of spin wave damping due to the spin-charge coupling mechanism as in metal- 
lic ferromagnets,— ^ which provides a quantitative understanding of the measured spin 
wave linewidth in neutron scattering studies of ferromagnetic manganites.— _ — The spin- 
charge coupling mechanism is especially important in saturated ferromagnets where the 
large inter-band gap forbids any low-energy spin wave damping due to decay into particle- 
hole excitations. Inelastic neutron scattering studies of iron pnictides indicate a constant 
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ratio r q /w q ~ 0.2 of the spin wave linewidth to energy, and clearly do not show any steep 
increase in the spin wave damping at higher energies as expected from decay into Stoner 
excitations.- 

In this paper, we will therefore investigate these two mechanisms for spin wave damping 
in the (0, 7r) ordered SDW state of itinerant electron models. We will also extend our earlier 
spin wave analysis in the (0, ir, tt) ordered SDW state^^ to include hopping anisotropy as 
appropriate for a layered system, and compare with the measured spin wave dispersion for 
iron pnictides. For this purpose, we will consider the single-band Hubbard model on a simple 
cubic lattice, with first and second neighbour hoppings t and t' in the x-y plane and a small 
interlayer hopping t z . 

We will also study the minimal two-band model with reference to the d xz and d yz orbitals 
of interest for iron pnictides,— with the effective hopping parameters ti - t 4 resulting from 
the hybridization of Fe 3d orbitals with themselves as well as through the As 3p orbitals 
lying above and below the square plaquettes formed by the Fe atoms. It is convenient to 
include the inter-orbital Hund's coupling term on an equal footing with the intra-orbital 
Hubbard interaction term within the general multi-orbital correlated electron model: 

H = — ^ ] tjj (O-i^o-O-jucr + a jva a i/J.v) ~ ^ ] U^Si^ • S{ u (1) 

where the interaction matrix elements U^ v = for /i = v and U^ v = 2J for fi ^ v refer to 
the intra-orbital and inter-orbital Coulomb interaction terms, respectively. For the two-band 
model, we will specifically consider the two cases: 

(i) ti=-l, t 2 =l, t3=t4=-0.3 at finite hole doping and 

(ii) ti=-l, £2=1-3, i3=t4=-0.85 at finite electron doping 

both of which yield stable metallic SDW states on including a finite Hund's coupling.— Case 
(ii) yields circular electron and hole pockets near the bottom of the upper band.— 

II. SPIN WAVE DISPERSION 

Figure 1 shows the calculated spin wave dispersion w q for the Hubbard model with 
anisotropic hopping, along with the neutron scattering experiment data points.- Here H,K,L 
correspond to wavevector q along the x,y,z directions. The overall energy scale is fitted with 
£=200 meV, and the interlayer-to-planar hopping ratio is taken as r=t z /t=0A. The other 
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FIG. 1: The calculated spin wave dispersion in the (7r,0,7r) ordered metallic SDW state of the 
Hubbard model with a small inter-layer hopping term and comparison with the neutron scattering 
data (from Ref. [4]). 

Hubbard model parameters are U/t w 11 (corresponding to A/t = 3 where 2 A = mil), 
t'/t = 0.3, and hole doping x ~ 40%, as earlier.— 

Analytical expressions for spin wave dispersion in the (0, 7r) and (0, 7r, 7t) ordered SDW 
states were derived earlier in two and three dimensions in the strong coupling limit of the 
t-t' Hubbard model at half filling.— It is useful to extend these expressions to approximately 
include the ferromagnetic coupling induced at finite doping, the finite sublattice magne- 
tization m in the doped state, and anisotropic hoppings in different directions as in iron 
pnictides. For the (0, 7r) state we thus obtain: 



2J r 
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, 2J ' 
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where J and J' represent first and second neighbour AF spin couplings in the metallic 
SDW state, which reduce to 4t 2 /U and 4t' 2 /U, respectively, in the insulating limit.— The 
coefficient b < represents the negative contribution of the carrier-induced ferromagnetic 
spin couplings in the x direction, whereas b = 1 in the insulating limit. Similarly, for the 
three-dimensional (tt, 0, ir) SDW state of a layered system, we obtain: 
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(3) 
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where r = t z /t is a small interlayer-to-planar hopping ratio. The energy scale on the left 
should be J(l + r 2 ) in general. The resulting spin-wave dispersions evaluated from Eqs. (2) 



FIG. 2: Spin wave dispersions w q (in units of J/m) evaluated from the analytical expressions given 
in Eqs. (2) and (3) for the (a) (0, tt) and (b) (ir, 0, n) ordered SDW states, obtained by including 
a doping-induced ferromagnetic spin coupling term and a small interlayer hopping term. 



and (3) are shown in Fig. 2 with 2J'/J = 2.0 and b = —0.3 for the two-dimensional case, 
and 2 J'/ J = 1.0, b = —0.3, and r 2 = 0.1 for the layered system. The dispersion in the latter 
case provides an excellent description of the neutron scattering results for iron pnictides. 

As seen from Fig. 2, the spin wave dispersions calculated from the approximate ex- 
pressions given in Eqs. (2,3) are remarkably close to that obtained from the t-t' Hubbard 
model, indicating that the essential features of the carrier-induced spin couplings are appro- 
priately incorporated. The substantial enhancement of J'/ J as compared to (t' /t) 2 reflects 
the additional second-neighbour AF spin couplings generated at finite doping. 

For the general spin wave propagator structure in the two-sublattice basis, including both 
retarded and advanced poles, we correspondingly obtain: 



x + (q,w) = -\ ( — 

2 V^q 



oj — w q + ii] u + ujq — ir] 



(4) 



a q -m(w/J) -/3 q 
— Pq a n + m(u/J) 

where the momentum-dependent terms a q and /3 q represent the terms in the two curly 
brackets in Eqs. (2) and (3). Including the magnetization m as above ensures that the spin 
commutation property ([S + , S~]) = (2S Z ) relating local spin correlations and magnetization 
is identically satisfied. 
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FIG. 3: The frequency dependence of the imaginary part of the particle-hole propagator x (qjw), 
including both intra-band and inter-band particle-hole excitations. The imaginary part vanishes 
over a large frequency range which includes the spin wave energy for mode q, implying vanishing 
contribution to spin wave damping due to decay into particle-hole excitations. 

III. SPIN- WAVE DAMPING DUE TO DECAY INTO INTRA-BAND 
PARTICLE-HOLE EXCITATIONS 

As mentioned earlier, spin-wave decay into low-energy intra-band particle-hole excitations 
in principle provide a mechansim for spin wave damping in a metallic SDW state. However, 
whether the simultaneous requirements of momentum and energy conservation afford suffi- 
cient phase space to yield significant linewidth is a relevant question. In the following, we 
will investigate this aspect from the imaginary part of for the t-t' Hubbard model, 

as well as for the minimal two-band models, both hole doped and electron doped. 

Fig. 3 shows the frequency dependence of the imaginary part of the particle-hole propa- 
gator u ) f° r the t-t' Hubbard model, including both intra-band and inter-band particle- 
hole excitations. In the two typical cases shown here, with wavevector q = tt/2 in the F and 
AF directions, the imaginary part vanishes over a large frequency range, which includes the 
spin wave energies for both modes (shown by arrows), implying vanishing contribution to 
spin wave damping due to decay into particle-hole excitations. Similar behaviour is obtained 
for the two-band model with £3 = £4 = —0.3, again indicating no intra-band contribution to 
spin wave linewidth. 

For the two-band model with £3 = £4 = —0.85, which yields circular electron and hole 
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FIG. 4: The spin-wave spectral function for the two-band model with circular electron and hole 
pockets,— showing significant broadening for wavevector q in the AF direction due to spin-wave 
decay into intra-band particle- hole excitations. 

pockets,— there is again no finite damping for q in the F direction. However, for q in the 
AF direction, we find substantial imaginary part at uj = w q , indicating available phase space 
for spin wave decay into intra-band particle-hole excitations. Fig. 4 shows the spin-wave 
spectral function in the (0, tt) state for three different wavevectors. The spin wave linewidth 
T q (half-width-at-half-maximum) is seen to increase with spin wave energy, and we obtain 
r q /w q ~ 0.25. Here the spin-wave spectral function was obtained from the RPA-level 
propagator as: 

from the complex eigenvalue A q (co>) (corresponding to the spin wave mode) of the matrix 
[C/][x°(q, a;)] involving the interaction matrix [U] given below Eq. (1) and the particle-hole 
propagator [x°(q, u)] in the two-orbital, two-sublattice basis.— 



IV. SPIN-CHARGE COUPLING MECHANISM FOR SPIN- WAVE DAMPING 

In view of the vanishing intra-band contribution to spin-wave linewidth obtained above 
for wavevector in the F direction due to absence of phase space for decay into particle-hole 
excitations, we will next investigate the contribution to spin-wave damping from the spin- 
charge coupling mechanism which is important in metallic ferromagnets. Within a system- 
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atic expansion of the irreducible particle-hole propagator beyond the RPA, the first order 
quantum corrections involving self-energy and vertex corrections due to electron-magnon 
coupling physically correspond to a second-order Raman scattering process involving spin- 
wave decay into a longer-wavelength intermediate-state spin wave accompanied with internal 
(majority-spin) charge excitations.— 1 ^ Effects of this spin-charge coupling mechanism both 
on the spin wave energy renormalization and spin wave damping have been studied in de- 
tail, and the calculated result r q /u; q ~ 0.1 is in good agreement with the inelastic neutron 
scattering measurements of spin wave excitations in ferromagnetic manganites.— ~— 

The first-order quantum correction to the irreducible particle-hole propagator (AA term 
in the sublattice basis) due to the spin-charge coupling mechanism is obtained as: 

0« (q, w) = J2 ™Q F sp-ch g p t-\ ~ — • ( 6 ) 

in terms of the SDW-state electronic energies The spin-charge coupling vertex: 

1 1 



sp— chg 



u 



Eg q - El e + u - i V Ef - El e + u - i V 



(7) 



identically vanishes at q = (0, 0) and (0, 7r) to ensure that the Goldstone mode is explicitly 
preserved. Here the A-sublattice magnon amplitude: 

1 / Jd!Q 



follows from Eq. (4) in terms of the magnon energy uq in the intermediate state. As 
the magnon spectrum is dominated by zone boundary modes with uq ~ Jo:q/to, this 
amplitude approximately reduces to the magnetization m as in metallic ferromagnets. Also, 
the majority amplitudes in the SDW state are approximated as 1, and only the corresponding 
dominant terms — interband terms (indicated by © and 0) in the interaction vertex and 
intra-band particle-hole terms (indicated by + and -) in the charge excitation propagator 
— were included in the quantum correction. 

The corresponding spin-wave linewidth due to spin-charge coupling is then obtained from 
the imaginary part of the self energy S sw (q, u) = mL r2 ( - 1 - ) (q, to), evaluated at the spin wave 
energy u = — w q corresponding to the advanced pole. The g-dependence of the spin-wave 
linewidth T q is shown in Fig. [5] for the (0, it) ordered SDW state of the t-t' Hubbard model 
(with same parameters A/t = 3, t'/t = 0.3, x ~ 40% as earlier). Except in the AF direction 
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FIG. 5: Momentum dependence of the spin-wave linewidth due to the spin-charge coupling mech- 
anism. Except for the AF direction starting from the zone center, the linewidth is around 10-20% 
of the spin wave energy over most of the Brillouin zone. 

from the zone center, over most of the Brillouin zone, the calculated linewidth is about 10- 
20% of the spin wave energy, in qualitative agreement with neutron scattering studies of iron 
pnictides. 4 Thus, while spin wave damping in the AF direction is dominated by intra-band 
particle-hole excitations, it is the spin-charge coupling mechanism which is important for 
modes in the F direction. 



V. CONCLUSIONS 



In conclusion, spin wave dispersion and damping in iron pnictides can be well understood 
in terms of physical mechanisms characteristic of metallic magnets, as shown by our com- 
parison with inelastic neutron scattering experiments. The carrier-induced ferromagnetic 
spins couplings as in metallic ferromagnets are crucial in stabilizing the (0, tt) and (0, tt, tt) 
ordered metallic SDW states, and analytical expressions for the spin wave dispersion were 
obtained in two and three dimensions by incorporating these spin couplings in expressions 
available for the insulating SDW states. On including a small interlayer hopping term, the 
calculated spin wave dispersion in the (tt, 0, tt) ordered SDW state was obtained in excellent 
agreement with experiments. 

With regard to spin wave damping, the intra-band particle-hole excitations and the spin- 
charge coupling mechanism presented distinctly complementary importance for modes in 
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the AF and F directions in the Brillouin zone. While for modes in the AF direction, the 
low-energy intra-band particle-hole excitations were found to yield significant linewidth for 
the two-band model with r q /a;q ~ 1/4, negligible contribution to linewidth was obtained 
for modes in the F direction, owing to severe phase-space restriction set by the simultaneous 
requirements of energy-momentum conservation. On the other hand, the spin-charge cou- 
pling mechanism, which is also important in ferromagnetic manganites, was found to yield 
significant spin wave damping for modes in the F direction, with Tq/cUq ~ 10 — 20%, in 
qualitative agreement with experiments. 

Evidence of Fermi surface folding associated with the SDW state has been observed in 
recent ARPES studies. 2 - Electronic quasiparticle dispersion and spectral function renor- 
malization in the SDW states due to electron-magnon interaction and multiple magnon 
emission-absorption processes should therefore be of interest. 
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